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The climate models that include a carbon-cycle need the vertical diffusivity of a passive tracer. Since an 
expression for the latter is not available, it has been common practice to identify it with that of salt. The 
identification is questionable since T, S are active, not passive tracers. We present the first derivation of 
the diffusivity of a passive tracer in terms of Ri (Richardson number) and R p (density ratio, ratio of salinity 
over temperature z-gradients). The following results have emerged: 
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(a) The passive tracer diffusivity is an algebraic function of Ri, R p . 

(b) In doubly stable regimes (DS, 0T/0z > 0 , 0S/0z< 0 ), the passive scalar diffusivity is nearly the same 
as that of salt/heat for any values of R p < 0 and Ri > 0 . 

(c) In DC regimes (diffusive convection, 07/0z < 0 , 0S/0Z < 0 , R p > 1 ), the passive scalar diffusivity is lar- 
ger than that of salt. At Ri = 0(1 ), it can be more than twice as large. 

(d) In SF regimes (salt fingers, 0T/0Z > 0 , 0S/0Z > 0 , R p < 1 ), the passive scalar diffusivity is smaller than 
that of salt. At Ri = 0(1 ), it can be less than half of it. 

(e) The passive tracer diffusivity predicted at the location of NATRE (North Atlantic Tracer Release 
Experiment) is discussed. 

(f) Perhaps the most relevant conclusion is that the common identification of the tracer diffusivity 
with that of salt is valid only in DS regimes. In the Southern Ocean, where there is the largest 
C0 2 absorption, the dominant regime is diffusive convection discussed in (c) above. 

Published by Elsevier Ltd. 


1. Formulation of the problem 

Since ocean general circulation models (OGCMs) cannot resolve 
small scale turbulent processes, the latter must be parameterized 
in terms of the resolved fields. OGCMs solve the time dependent 
equations for the mean velocity, mean temperature, mean salinity 
and mean concentration of passive tracers such as C0 2 , since the 
carbon cycle is an important diagnostic in climate studies. Some 
climate studies also include CFC. While both temperature and 
salinity are active tracers, for they affect the velocity field via the 
buoyancy term and thus alter the status of the turbulent flow, trac- 
ers such as C0 2 , CFC, Oxygen, if sufficiently diluted, as is ordinarily 
the case in the ocean, are passively carried along by the flow with- 
out altering it. * 1 It follows that the diffusivity of a passive tracer must 


* Corresponding author at: NASA, Goddard Institute for Space Studies, New York, 
NY 10025, USA. Tel.: +1 212 6785571; fax: +1 212 6785560. 

E-mail address: vcanuto@giss.nasa.gov (V.M. Canuto). 

1 As pointed out by a referee, there are cases of lakes in which the concentrations of 
C0 2 and CH 4 contribute significantly to the density stratification to give rise to 
double-diffusive convection (Schmid et al., 2010). 
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be different from that of heat and salt. Identifying the passive tracer 
diffusivity with that of salt is, however, a common practice probably 
because no expression for a passive tracer diffusivity is available in 
the literature. 

In a recent work (Canuto et al., 2010, CIO), 2 we derived the ver- 
tical diffusivities for momentum, heat and salt as functions of the 
two large scale parameters Ri and R p (Richardson number and the 
density ratio) defined as follows: 


where £ = (2SySy) 1/2 is the mean shear with 2Sy = U,j + U j}i with i, 
j= 1, 2, 3; T,S and U are the mean potential temperature, mean 
salinity and mean velocity fields; we employ the notation a z = da / 
0z. The thermal expansion and haline contraction coefficients 
oc,P = (-p _1 0p/0T, +p _1 0p/0S) are computed using the non-linear 
UNESCO equation of state (UNESCO, 1981). The density ratio R p 


2 In order to make this presentation self-contained, in Appendix A we report the 
expressions derived in CIO which are needed here. 
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accounts for double-diffusion (DD) processes such as SF -salt fingers , 
0 < R p < 1 (warm-salty water over cold-fresh water), DC-diffusive 
convection, R p > 1 (cold-fresh water over warm-salty water) and 
DS -doubly stable R p < 0 (warm-fresh water over cold-salty water). 
Specifically, in CIO it was shown that the solutions of the Reynolds 
stress model (RSM) yield the following expressions for the 
diffusivities where the subscript a stands for momentum, heat 
and salt: 

i< x = ^fs x = r*, r x = }rS. x (xN ) 2 (2) 

Here, I< is the eddy kinetic energy, e is its rate of dissipation, 
t = 2 Kir 1 is the dynamical time scale and N is the Brunt-Vaisala 
frequency, N 2 = -gp^dp/dz. The functions r a are known as mixing 
efficiencies and S a are dimensionless structure functions both of 
which are given by the RSM as functions of: 

r x [Ri,R p ,(xN) 2 }, S a [Ri,R p ,( tN) 2 ] (3) 

In the mixed layer the first representation in (2) is more conve- 
nient since one can use a I<-s model to determine these two vari- 
ables while in the thermocline the £-equation is quite uncertain 
and the second representation in (2) is more convenient. In fact, 
IGW (internal gravity waves) studies (Polzin et al., 1995; Polzin, 
1996; Kunze and Sanford, 1996; Gregg et al., 1996; Toole, 1998), 
provide an expression for s N~ 2 . The task of the turbulence models 
is that of modeling the mixing efficiencies r a ( Ri , R p ). 

Since in CIO we did not derive an expression for the diffusivity 
of a passive tracer, here we remedy that shortcoming by presenting 
a derivation for such a diffusivity. 

The derivation is not trivial not because of mathematical but 
rather physical reasons and that is why we include it in the main 
text rather than in an Appendix. The difficulties are of two kinds. 
First, since the starting Eq. (4) is valid for an arbitrary tracer, during 
the derivation one must identify where and how a generic tracer is 
treated as passive or active. This is a somewhat delicate problem 
but not a very difficult one. The most difficult part is the determina- 
tion of the dissipation time scales that characterize small scale pro- 
cesses in the passive tracer case. In CIO, we presented a lengthy 
discussion of the analogous difficulty when one treats active trac- 
ers such as temperature and salinity. Assuming a constant ratio of 
the dissipation time scales to the dynamical time scale, as it was 
assumed for many years, ultimately because lack of data did not 
motivate other suggestions, yielded a rather poor fit to the NATRE 
data (North Atlantic Tracer Release Experiment, Ledwell et al., 
1993, 1998; St. Laurent and Schmitt, 1999, see Fig. 3 of Canuto 
et al., 2008). A better fit was achieved when the time scale ratios 
were made Ri, R p dependent. Within the present context, the ques- 
tion is how to relate those dissipation times scales to the new ones 
that characterize a passive tracer. As we discuss below, the second 
part of the derivation consists of finding such a relation which we 
believe we have done without introducing arbitrary or heuristic 
assumptions. Rather, we required that the diffusivity that we ob- 
tain with a direct solution of the turbulent equation for a passive 
scalar, ought to coincide with those of either salinity and/or tem- 
perature when the latter are considered so diluted as not to affect 
the density field and thus behave like passive tracers. Such a 
requirement allows to find a relation between the dissipation time 
scale of a passive tracer and those already determined in the tem- 
perature-salinity case. 

2. Passive tracer vertical flux J c = w c' 

We begin with the general equation for an arbitrary tracer (of- 

ten referred to as a scalar) which at this point can be either passive 

or active. Such an equation is 


dtpc + di(pcui) = x c pS?c, (4) 

where the rhs represents molecular effects with a kinematic diffu- 
sivity Xc- Following standard procedure, the total field c is split into 
an average and a fluctuating component, c = C + c! and the same is 
done with the total velocity field U; = U, + u\. Next, we introduce 
the above decompositions in (4), use the incompressibility condi- 
tion, average the result and subtract it from the equation for the total 
c-field. These purely algebraic steps lead to the following equation 
for the fluctuating component of the generic field c! ( D t = S t + U^i): 

Dtd + u'fiiC + a,- ( u'id - up) = Xctfc' (5) 

Using an analogous procedure for the velocity field described by the 
Navier-Stokes equations, one derives the equation for the field w' 
which reads as follows: 

s t w' = -p 0 1 0 z p' - gp Q 1 p' - dj (w'Uj - wq) + vS ifW (6a) 

where v is the kinematic viscosity of seawater. As expected, the 
averages of both (5)— (6) yield the identity 0 = 0 since all fluctuating 
variables have zero average. 

Next, we recall that if the field c is not affecting the density, and 
this is the assumption about c being a passive tracer, then the den- 
sity is only affected by the truly active scalars, temperature and 
salinity and thus the density field p' in (6a) is taken to be: 

p'Po 1 = -cuff + fis' (6b) 

where O', s' are the fluctuating components of the T and S fields. The 
next step is again purely algebraic and consists of multiplying Eq. 
(5) by w' and Eq. (6a) by c', averaging the results and summing 
the two equations. The result is the dynamic equation for the verti- 
cal flux of the passive tracer J c = Wc': 

0 J c = - w^0 z C + g (u0e - jSS 7 ?) - r pc } J c (7) 

where we have lumped together the pressure correlation and the 
dissipation terms in the dissipation-relaxation time scale x pc . Next, 
we must find an expression for the two covariances 0 / c / ,s 7 c 7 that en- 
ter Eq. (7). This is done by using the two equations for the fluctuat- 
ing temperature and salinity fields which are obtained using a 
procedure analogous to the one used to arrive at (5). The results, 
which have been derived in the literature in many textbooks (e.g., 
Leslie, 1973; Lesieur, 1990), are: 

D t e' + UjdjT + 0,- (u;e' - ufi) = x^fd' (8a) 

D t s' + U'fiiS + di ( U-S ' - ujs 7 ) = Xs d f S ' (8b) 

where Xe,s denote the heat and salt molecular diffusivities. We re- 
call that Xe ~ 100Xs and that the ratio v/xo « 7 is the Prandtl num- 
ber for seawater. Multiply (8a) by c' and (5) by O', sum the two 
relations after averaging them and carry out an analogous operation 
with (8b). The result is the two equations: 


0 t 0'c' = - 1 

{jhC,z +JJ*) - 'C C dO' c ' 

(9a) 

dtWc 7 = - ( 

JsQ +J c s,z) - T “W 

(9b) 


where, as in Eq. (7), the two time scales t c9 and t cs represent the 
relaxation-dissipation time scales of the corresponding correlations 
(for simplicity, we have neglected the third-order moments). 

Before we proceed with the derivation, we pause momentarily 
to point out an important difference between passive and active 
tracers. Notice that Eqs. (9a) and (9b) are independent of the den- 
sity and they are therefore quite general. Suppose however that c! 
had been an active scalar: if c! had been either O', s', Eqs. (9a) and 
(9b) would have given rise to the terms: 
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O' 2 , S' 2 (9c) 

as expected since they represent the potential energies stored in the 
0, s turbulent fields which are part of the energy budget; however, 
there is no such a term as c' 2 entering the above equations (even 
though it may be constructed from Eq. (5)) since a passive scalar 
field does not contribute to (affect) the energy balance. In the der- 
ivation of the mixing model presented in CIO, the dynamic equa- 
tions for (9c) were derived in Appendix A, Eqs. (A3) and (A4) of 
CIO. In the c' -► O' case, for example, Eq. (7) would become the 
equation for the heat flux and the second term would be twice 
the potential energy in the temperature field, representing a source 
of the heat flux. In such a case, Eq. (7) would coincide with Eq. (Al) 
of CIO. 

In Eqs. (9a) and (9b), Jh, s are the vertical fluxes of heat and salt 
whose form was derived in CIO to be: 


k = WO' = - I< h T , z , J s = Ws' = -I( s S z (9d) 

where the heat and salt diffusivities K hs are of the general form (2) 
and the corresponding r h , s are given in Appendix A. Taking the sta- 
tionary limits in Eqs. (7), (9a), (9b) and using (9d), we finally obtain 
the desired form of the vertical flux of a passive tracer: 


j C = w'c = -k c ^, K c = r c 3- ) 
Jc C 0Z’ C N 2 


,, (10a) 

where the mixing efficiency for the passive tracer is derived to be: 


r c = 


Pi(tN) 2 


1 +P 1 P 2 (tN) 


w 




2 \ 2K 


(10b) 


Relation (10b) deserves several comments. First, r p is the mix- 
ing efficiency of the buoyancy (density) field and is defined in 
terms of the heat-salt mixing efficiencies J\ s via the following 
relation: 


r p = (r h -R p r s )(\-R P r\ 


(10c) 


Physically, (13b) means that salinity is so diluted as not to affect the 
density, in which case salinity behaves like a passive scalar. In such 
a limit, we therefore expect that: 

R p ^0: I( c — > I( s . (13c) 

Furthermore, from (A5) and (A9) we have the following relations: 
Rp ^ 0: r^oo, P — > 7l 4 715 - 7T47T2 , q — > 7Ti 7T 2 - 

(13d) 


At the same time, since the function A p has the same structure as 
(10c) w.r.t. heat and salt, we have the two additional relations: 

x-W'O-V <13e) 


where we used the first of (A4) and (13d). Introducing (13e) into 
(13a), a short algebra leads to the final form: 


n A . p 1 [l+7T 4 (7r 5 -p 2 )x] 

p ' c (1 +7r 4 7r 5 x)(l +P,p 2 x) 


(13f) 


Next, we show how (13f) is related to the salt case. To do so, we em- 
ploy both relations (A4) and (A5) and obtain: 

a = Ah = A h = 7Ti (1 + px) 
s rR p K h /K s h n 4 (1 +qx) 

71) [1 +7r 4 (7r 5 -n 2 )x] 

(1 + 7T 4 7r 5 X)(l + 7Ti tl 2 x) ’ ^ ^ 

where in the last step we used (13d) forp and q. Eq. (13g) coincides 
with (13f) provided we identify: 

Rp — > 0 : p, =7t,(Rp = 0), p 2 =n 2 (Rp=0). (13h) 

Next, we treat the case when heat becomes passive which is repre- 
sented by the following relations: 

a -> 0, R p —> oo : I( c -* I< h (14a) 

In this limit, we have the following relations: 


where the r htS are given by Eqs. (A1)-(A5) and (A9). Second, the 
determination of the variable (t N) 2 was discussed at length in Sec- 
tion 4.4 of CIO (for the P = s case) and we need only mention the fi- 
nal result presented in Eqs. (A10)-(A14). Third, the functions pi i2 
are defined as follows: 




T cs 

T’ 


al) 


whose determination is the most difficult part of the problem as 
discussed below. We recall that in CIO, we introduced the following 
dimensionless time scales: 





(12a) 


which are given in terms of Ri and R p by relations (A13), (A14) 
which were assessed against NATRE data in Canuto et al. (2008, 
C8). Since relations (A13) show that the forms of conserve the 
heat-salt symmetry, there is no a priori justification to identify p ? 
with either since in principle both are potentially acceptable 
candidates. The construction of the function: 


Pi =/(7Ti, 7T 4 ) 


(12b) 


P — * —712^4 Rp, q — * 7Ti (7T 2 — n 3 )Rp, 


x->~x 2 N 2 R~\ 


px -> n 2 n 4 T 2 N 2 , qx -n\ (n 2 - n 3 )x i N i , 


2 ^r2 


^ _ I< h _ n 4 1 + qx 7T 4 1 — (n 2 - n 3 )z 2 N 2 

P ~K s ~TX l TTpx^ Th 1 + 7t 2 7t 4 T 2 N 2 


The function A h defined in (A4) becomes: 
7l 4 


Ah 


1 +7l 2 7l 4 1 + 


K)' 1 ] 


T 2 N 2 


(14b) 

(14c) 


or, using the last of (14b) and after several algebraic steps, we 
obtain: 


Ah 


7l 4 [l - 7Ti (7l 2 - 7T 3 )T 2 N 2 j 


1 + (7Ti 7T 3 + n 2 n 4 )T 2 N 2 + 7ri7r 2 7i 3 7i 4 T 4 N 4 

On the other hand, we also have that: 


A P 


Ah_ 

rRp ' 


(14d) 


(14e) 


is a problem we tackle next. To do so, we use (Al) to transform 
(10b) to the following form: 




Pi 

1 +p 1 p 2 (TN) 2 


1-p 2 (TN) 2 Ap\. 


(13a) 


Next, consider the active tracer salt in the limit in which the haline 
contraction coefficient becomes vanishingly small: 

P - (0p/0S ) T ^0, R p -> 0. (13b) 


Using (14d) and (14e), we now have that (13a) becomes: 

Pi [l - n 3 (p 2 - u 3 )t 2 N 2 ] 

J\ c L J (14f) 

1 + (n 3 n 3 +p,p 2 )t 2 N 2 + n^n 3 p^p 2 -c 4 N 4 

Relations (14f) and (14d) coincide provided we identify: 

a^O, Rp^ oo: p 1 = 7l 4 (Rp = oo), p 2 = 7l 2 (Rp = oo), (14g) 

which is the second limit we shall use to construct relation (12b). 
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Using the two relations (13h) and (14g), we suggest the follow- 
ing form for (12b): 

Rp + ft/ 


p 2 = n 2 . 


(17) 


Pi =/(7C 1 ,7C 4 ) = - 


(15) 


+Rp 

which we arrived at considering that since the shortest time scale 
dominates, one begins with the relation: 

J_ _ A_ B_ 

T nc 


(16) 


L pc L ps T'pO 

Using the tt’s from (12a) and pi from (11), we chose A, B to satisfy 
the two requirements (13h) and (14g). The result is (15). As for 
p 2 , its identification with n 2 seems unproblematic: 


Using the variables (tN) 2 , r p , pi >2 and (A2), (A3), in Figs. 1-3 we 
plot r c given by (10b) for different values of Ri and R p . 

3 . Results 

In Fig. 1, corresponding to the doubly-stable, DS case, R p < 0, we 
present the dimensionless mixing efficiencies for a passive scalar, 
panel (a) and for spice, panel (b). The latter is defined as a = a 
6' + p s' and since it is invariant under the exchange of the T and 
S fields, it might be considered as a possible candidate for a passive 
tracer. The corresponding mixing efficiency is given by 


DS: R p =0(. 
(a) r c 


— 0 . 1 (— ), - 1 (_), — 10 (— 

(b) r CT 
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(d) r c /r h 
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Fig. 1. Doubly stable case, 5T/6z> O,0S/0z< 0, R p < 0. The dimensionless mixing efficiencies of a passive tracer and spice and the ratios of the former with salt, heat, 
buoyancy and spice. Spice is defined as o = o tO 1 + /3s 
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DC: R„=oo( ), 5( ), i.67( ), 1.25(-.-.) 


(a) r c 


(b) r a 







Fig. 2. Same as in Fig. 1 but for the case of diffusive convection 5 T /dz < 0, dS/dz < 0, R p > 1. The values of R p are just the reciprocal of those of the SF case shown in Fig. 3. 


r a = (r h + R p r s )(\+R p y' ( 18 ) 

In the remaining panels (c)-(f) we exhibit the ratios of the passive 
tracer with salt, heat, buoyancy and spice. As one can observe from 
panels (c)-(f), to all practical purposes, there is little difference be- 
tween the diffusivity of a passive scalar and that of salt, heat or 
spice. In principle, this might be construed as a justification of the 
common assumption that a passive tracer diffusivity can be identi- 
fied with that of salt. However, the Southern Ocean which is the 
largest C0 2 sink due to its largest extent (Takahashi et al., 2009), 
is not primarily characterized by R p < 0 but by diffusive convection, 
DC, R p > 1. Therefore, in carbon-cycle studies, identifying the pas- 
sive tracer with salt may lead to incorrect results. The DC results 


are presented in Fig. 2. As one can observe from panel (c), the pas- 
sive tracer diffusivity is larger than that of salt up to a factor of two 
and it is smaller than that of heat by the same amount, panel (d). 
Finally, in Fig. 3 we present the results for the SF case, in which case 
we notice that the tracer mixing efficiency is smaller than that of 
salt but larger than that of heat by roughly the same amount. 

4. Assessment of the model prediction against the NATRE data 

Ledwell et al. (1993, 1998) reported passive tracer diffusivity 
measurements during the NATRE (26 N, 28 W) experiment. In 
Fig. 4 we reproduce their data (open triangles) at a depth of 
300 m for two successive periods (the point on the left corresponds 
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SF: R p = 0( ), 0.2( ), 0.6( ), 0.8(-.-.) 


(a) r c 


(b) r CT 




Fig. 3. Same as in Fig. 1 salt fingers, dT/dz > 0, dS/dz > 0, R p < 1. 


to the first 6 months’ measurements while the one on the right cor- 
responds to the subsequent 24 months). The model results (de- 
noted by a cross) are closer to the data of the first measurement. 
This may be due to the fact that the values of R p used to compute 
the model results (see description below) were obtained from St. 
Laurent and Schmitt (1999), cited as SS99 and corresponded to a 
time closer to the measurement of the first period. Since we have 
no data other than those at 300 m, the z-profiles must be viewed 
as theoretical predictions to be assessed when future measure- 
ments will become available. The reason why the model results 
(crosses) exhibit an error bar is because they are based on the mea- 
sured values of R p which have errors. 

As suggested by a referee, in the same figure we also plot the z- 
profile of the salt diffusivity I< s (black squares). As in the tracer 


case, the model salt diffusivities have an error bar because they de- 
pend on the measurements based values of R p . As expected on the 
grounds of Fig. 3(c), in a salt-fingers regime such a one at NATRE, 
one expects to have I< s > I< c which is indeed the case in Fig. 4. 

It is important to discuss how the model results (crosses and 
black squares) were obtained without the use of an OGCM. We em- 
ployed the following procedure. The z-profile of the density ratio 
R p was taken from the NATRE data from SS99. The functions A htS 
were constructed using (A4) in which the heat to salt flux ratio 
“r” was computed from its definition (A5). The variable x, defined 
in (A9), was obtained by solving Eqs. (A10), (All), (A12); the var- 
iable w' 2 /I< is given by Eqs. (A.2) and (A3). As for the Richardson 
number, while in CIO its background value Ri(bg) in the thermo- 
cline was taken to be 0.5, we have since then found that an R p 
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Fig. 4. The z-profile of the model passive tracer and salinity diffusivities and the NATRE measurements. For details, see Section 4. 


x passive tracer, model 
■ salt, model 

a passive tracer, measurement 
i— i error bars 


dependent Ri{bg) yields a better fit to the SS99 data on heat and 
salt diffusivities. Since we believe that Ri{bg) is set by a self-limit- 
ing process in which the turbulent viscosity produced by the wave 
breaking interferes with the waves themselves, we tried the 
hypothesis that in the background r m , which is proportional to 
the viscosity, is a constant. Since r m is a function of Ri and R p , 
adopting this hypothesis yields Ri{bg) vs. R p . The best fit to the 
SS99 data was obtained with r m = 0.46 and the Ri(bg) correspond- 
ing to that value was accordingly used. 


5. Maps of active and passive tracer diffusivities 

In Figs. 5-10 we present meridional maps, depth vs. latitude (at 
a longitude of 29 W passing through the NATRE location), of the 
vertical diffusivities (cm 2 s _1 ) for the following fields: 

momentum, heat, salt, spice, buoyancy, passive scalar. 
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Fig. 5. Depth-latitude maps of the momentum diffusivity at 29 W (near the NATRE 
location). For comments see Section 5. 


Fig. 6. Same as Fig. 5 but for the heat diffusivity. 

GISS 



Fig. 7. Same as Fig. 5 but for the salt diffusivity. 
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GISS 



Fig. 8. Same as Fig. 5 but for the spice o = oc 6' + /3 s' diffusivity. 


GISS 



Fig. 10. Same as Fig. 5 but for the passive scalar diffusivity. 


GISS 



Fig. 9. Same as Fig. 5 but for the buoyancy diffusivity. 
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Fig. 11. Ratio of the passive scalar to salt diffusivity in the same region as in Figs. 5- 
10. For comments see Section 5. 


The vertical mixing model is given in CIO except that the Ri{bg) was 
determined as discussed in Section 4. The inputs to the mixing 
model (velocity, temperature and salinity fields) were obtained 
from a snapshot of existing model output from a coupled ocean- 
atmosphere general circulation model, specifically, the GISS coupled 
Model E with 1° x 1.25° horizontal resolution and 32 vertical levels. 
The snapshot was taken after 150 years of model spin-up in which 
the KPP vertical mixing scheme had been used. Several comments 
are in order: in all cases, the large values found near the top corre- 
spond to the enhanced diffusivities in the mixed layer. For all the 
diffusivities, the low values at the equator are due to the latitude 
dependence of the breaking of the IGW (internal gravity waves). 
As expected in a stably stratified case, the momentum diffusivity, 
Fig. 5, is larger than both heat Fig. 6 and salt Fig. 7 by a sizeable 
but not uniform, factor: the momentum diffusivities are between 
0.1 -0.3 while the heat diffusivity is at most 0.1. Salt diffusivities 
are larger than heat in places where SF (salt fingers) operate, as seen 
by the green patches in Fig. 7. Much of the effect of DD is however 


washed out in the case of spice, Fig. 8, whose diffusivity is more 
similar to that of heat than salt. In Fig. 9 we exhibit the diffusivity 
I< p , Eqs. (2) and (10c), corresponding to buoyancy. Since in regions 
of strong SF or DC such a diffusivity is negative, one can observe 
the patches where that occurs. For the vast majority of ocean re- 
gimes, the values are, as expected, smaller than those of heat and 
salt separately. In Fig. 10 we present the passive scalar diffusivity 
that can be compared with the heat and salt in Figs. 6 and 7. To 
make the comparison easier, in Figs. 11 and 12 we exhibit the ratios 
of the passive tracer diffusivity to those of salt and heat. The green 
regions in Fig. 11 ,k c <k s , and the red regions, k c > K h , in Fig. 12, cor- 
respond to a SF regime (see Eq. (19)). The yellow regions correspond 
to doubly stable regimes in which the passive scalar is very close to 
heat and salt, as seen in Fig. l(c)-(d). To further visualize the results 
of the model, in Figs. 13 and 14 we exhibit the momentum and heat 
diffusivities of the KPP model (Large et al., 1994) which can be com- 
pared with Figs. 5 and 6, respectively. The different texture of the 
results is apparent. 
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Fig. 12. Same as Fig. 11 but for the ratio of passive scalar to heat diffusivity. 



Fig. 14. Heat diffusivity in the KPP model. The different texture of the results from 
those in Fig. 6 is apparent. 
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as small as that of salt and twice as large as that of heat, a situation 
that can be summarized as follows: 

DS: K c ^I< h ^K s , 

DC: /C c >/C s , K c <K h , (19) 

SF: /C c </C s , K c >K h . 

The next step is to employ relation (10b), together with the vertical 
mixing scheme in CIO, in an OGCM with a carbon-cycle model and 
assess the difference with the previous assumption of a tracer rep- 
resented by salt. 
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Appendix A. Ingredients of the model 

General relations: 


Fig. 13. Momentum diffusivity in the KPP model. Notice the different texture of the 
results from Fig. 5. 


6. Conclusion 


1 9 w /2 W' 

r x = -(TN) 2 A x — , S a = A X — . (A1) 

Ratio w /2 /2K. From Eqs. (37) and (39) of CIO it follows that: 
General case: 


In this paper we have derived a general expression for the dif- 
fusivity of a passive tracer in terms of large scale variables Ri and 
R p . The result is given by relations (10a), (10b). The most relevant 
conclusion of this work is that the common identification of a tra- 
cer diffusivity with that of salt is valid only in regions of doubly 
stable (DS) T-S gradients which does not occur in much of the 
Southern Ocean where the largest C0 2 absorption occurs. In those 
regions, the dominant feature is an unstable T-gradient and a sta- 
ble S-gradient, a regime known as Diffusive Convection, DC, as the 
maps (courtesy of D.E. Kelly) shown in Figs. 6 and 7 of CIO show. 
The results in Fig. 2 show that in the DC regime the tracer diffusiv- 
ity for values of Ri of order unity, is roughly twice as large as that of 
salt and twice as small as that of heat. Conversely, in the SF case, 
the tracer diffusivity for values of Ri of order unity, is roughly twice 


w' z 

2 K 


\+^A p (TN) 2 +±A m (TZ) 2 


which, for the production = dissipation case, simplifies to: 


w 

P = £: 2 K = 


^-+\m 2 


Fleat and salt: 

Ah =7Lt[l +px + n 4 n 2 x(l -r” 1 )] -1 , A s =A h (rR p )~ 1 

Fleat-to-salt flux ratio: 


oc t w6 _ 1 I< h _n 4 1 1 + qx 
a s Ws ~ R p K s ~ 7Ti R p 1 + px ’ 


(A2) 


(A3) 


(A4) 

(A5) 
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Momentum: 


where a = 10. For R p < 0 and Ri > 0, we have: 


A m = 


Am 1 


Ami = 


714-71, + (ji 1 -^(l-r 1 ) 

1 


xAh 


A m2 = 10 + (7t 4 - n 1 R p )x + 

Dimensionless variables x, p and q (N 2 = -gp^dp/dz): 

x= (tN) 2 (1 -R p )~\ p = n 4 n 5 -n 4 n 2 (1 +R P ), q 
= 7T, 7T 2 (1 +R P ) - n^n 3 R p 
Cubic equation for G m = (tZ ) 2 = (rN) 2 Ri~ u . 


c 3 C m +c 2 C m 

with: 


■ Ci C m + 1—0, 


(A6) 

(A7) 

(A8) 

(A9) 

(A10) 


C 3 = A ] Ri i A 2 Ri 2 . c 2 = A 3 Ri 2 + A 4 Ri. c, = A 3 Ri + Ag , (All) 
where: 

150(1 - R p ) 3 Ai = nin 4 (n 4 - niR p ){n 2 (\5n 3 + 7) (r 2 p + l) 

+ [14(7r 2 -n 3 )-\5n 2 3 ]R p } 

9000(1 -R p ) 2 A 2 =Ui n 4 {n 2 {2\0ni - 1507i 3 + 7)(R 2 p + l) 

+ [14(7i 2 - 7 T 3 ) (1 + 15tti + 1 5tt 4 ) + 1 507r|]i? p + 2\0n 2 (n 4 - 7ii)} 
150(1 - R p ) 2 A 3 = 7ii[5n 2 n 4 (30n 3 + 17) + 7 Ti (15713 + 7)] 

(R 2 p + l) - (15713 + 7) (71? - 71^) - [10711713 714(15713 + 17) 

+ 157l 2 (7r? + 714 ) + 147T! 7T 4 (1 - 10tt 2 )]^ 

9000(1 - R p )A 4 = [1 50(7T! tt 3 + 7T 2 7r 4 ) - 77r 4 (1 + 307ii)] 

R p - 150(7ii7r 2 + 7T 3 7T 4 ) + 7tt 4 (1 +307i 4 ) 

30(1 - ^ p )A 5 == [—30(71! 7T 3 + n 2 n 4 ) - \7ni}R p 
+ 30(7Ti7i 2 + n 3 n 4 ) + 1 77T 4 , A 6 = -1/60 

(A12) 


7Ti )4 = 71? 4 (1 + Ri) \ 712,3,5 = 7r° 35 . (A14) 

For Ri < 0, 7i, < = 7i? for any k. In relations (A13) Ko = 1.6 is the Kol- 
mogorov constant. 
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The 7i’s: For R p ^ 0 and Ri ^ 0: 

tt 2 = 7I°(1 +Rir’ [l +2RiRp(\ +R 2 )“’ , 


Ri 

1 + ciR p J 


n 5 = n\, 


(A13) 


7i° = 7i° = (27Ko 3 /5) ' /2 (1 + 

7i? = 1 /3, n 3 = n° 3 = n° 5 = (j t = 0.72. 


